Mechanisms of long-term decline in size of lesser sandeels in the North Sea explored using a growth and phenology model

The lesser sandeel Ammodytes marinus is a key species in the North Sea ecosystem, transferring energy from planktonic producers to top predators. Previous studies have shown a long-term decline in the size of 0-group sandeels in the western North Sea, but they were unable to pinpoint the mechanism (later hatching, slower growth or changes in size-dependent mortality) or cause. To investigate the first 2 possibilities we combined 2 independent time series of sandeel size, namely data from chick-feeding Atlantic puffins Fratercula arctica and from the Continuous Plankton Recorder (CPR), in a novel statistical model implemented using Markov Chain Monte Carlo (MCMC). The model estimated annual mean length on 1 July, as well as hatching date and growth rate for sandeels from 1973 to 2006. Mean length-at-date declined by 22% over this period, corresponding to a 60% decrease in energy content, with a sharper decline since 2002. Up to the mid-1990s, the decline was associated with a trend towards later hatching. Subsequently, hatching became earlier again, and the continued trend towards smaller size appears to have been driven by lower growth rates, particularly in the most recent years, although we could not rule out changes in size-dependent mortality. Our findings point to major changes in key aspects of sandeel life history, which we consider are most likely due to direct and indirect temperature-related changes over a range of biotic factors, including the seasonal distribution of copepods and intra- and inter-specific competition with planktivorous fish. The results have implications both for the many predators of sandeels and for age and size of maturation in this aggregation of North Sea sandeels.

Patterns of jellyfish abundance in the North Atlantic

A number of explanations have been advanced to account for the increased frequency and intensity at which jellyfish (pelagic cnidarians and ctenophores) blooms are being observed, most of which have been locally directed. Here, we investigate seasonal and inter-annual patterns in abundance and distribution of jellyfish in the North Atlantic Ocean to determine if there have been any system-wide changes over the period 1946–2005, by analysing records of the presence of coelenterates from the Continuous Plankton Recorder (CPR) survey. Peaks in jellyfish abundance are strongly seasonal in both oceanic and shelf areas: oceanic populations have a mid-year peak that is more closely related to peaks in phyto- and zooplankton, whilst the later peak of shelf populations mirrors changes in SST and reflects processes of advection and aggregation. There have been large amplitude cycles in the abundance of oceanic and shelf jellyfish (although not synchronous) over the last 60 years, with a pronounced synchronous increase in abundance in both areas over the last 10 years. Inter-annual variations in jellyfish abundance in oceanic areas are related to zooplankton abundance and temperature changes, but not to the North Atlantic Oscillation or to a chlorophyll index. The long-term inter-annual abundance of jellyfish on the shelf could not be explained by any environmental variables investigated. As multi-decadal cycles and more recent increase in jellyfish were obvious in both oceanic and shelf areas, we conclude that these are likely to reflect an underlying climatic signal (and bottom-up control) rather than any change in fishing pressure (top-down control). Our results also highlight the role of the CPR data in investigating long-term changes in jellyfish, and suggest that the cnidarians sampled by the CPR are more likely to be holoplanktic hydrozoans and not the much larger meroplanktic scyphozoans as has been suggested previously.

Spatial and inter-decadal variability in plankton abundance and composition in the Northwest Atlantic (1958-2006)

The results of Continuous Plankton Recorder sampling in the NW Atlantic between 1958 and 2006 are presented for 11 plankton taxa in eight shelf and deep ocean regions. For shelf regions, phytoplankton abundances increased in the early 1990s, mainly in winter, as the contribution of Arctic-derived freshwater to the Newfoundland (NLS) and Scotian shelves (SS) increased. Farther east, in the sub-polar gyre, phytoplankton levels increased with rising temperatures during the 1990s and 2000s. In both areas, the changes can be explained by increased stratification. The increased influx of arctic water to the NLS in the 1990s was also probably directly responsible for the increased abundances of two arctic Calanus species (C. glacialis and C. hyperboreus) and indirectly responsible for the decreased abundance of Calanus I–IV (mainly C. finmarchicus), perhaps via changes in food composition. On the SS the arctic Calanus species increased in abundance in the 2000s, likely as the result of increased transport from the Arctic via the Gulf of St Lawrence. In the deep ocean, plankton seasonal cycles changed little over the decades and increasing phytoplankton levels in the 2000s were accompanied by increases in zooplankton abundance, suggesting bottom-up control. In shelf regions, phytoplankton increases in the 1990s were in winter and Calanus I–IV appeared earlier in spring than in previous decades. Zooplankton levels generally did not change overall however, perhaps because the species examined were mainly inactive during winter.

Environmental context explains Lévy and Brownian movement patterns of marine predators

An optimal search theory, the so-called Levy-flight foraging hypothesis(1), predicts that predators should adopt search strategies known as Levy flights where prey is sparse and distributed unpredictably, but that Brownian movement is sufficiently efficient for locating abundant prey(2-4). Empirical studies have generated controversy because the accuracy of statistical methods that have been used to identify Levy behaviour has recently been questioned(5,6). Consequently, whether foragers exhibit Levy flights in the wild remains unclear. Crucially, moreover, it has not been tested whether observed movement patterns across natural landscapes having different expected resource distributions conform to the theory's central predictions. Here we use maximum-likelihood methods to test for Levy patterns in relation to environmental gradients in the largest animal movement data set assembled for this purpose. Strong support was found for Levy search patterns across 14 species of open-ocean predatory fish (sharks, tuna, billfish and ocean sunfish), with some individuals switching between Levy and Brownian movement as they traversed different habitat types. We tested the spatial occurrence of these two principal patterns and found Levy behaviour to be associated with less productive waters (sparser prey) and Brownian movements to be associated with productive shelf or convergence-front habitats (abundant prey). These results are consistent with the Levy-flight foraging hypothesis(1,7), supporting the contention(8,9) that organism search strategies naturally evolved in such a way that they exploit optimal Levy patterns.

Modelling the effects of climate change on the distribution and production of marine fishes: accounting for trophic interactions in a dynamic bioclimate envelope model

Climate change has already altered the distribution of marine fishes. Future predictions of fish distributions and catches based on bioclimate envelope models are available, but to date they have not considered interspecific interactions. We address this by combining the species-based Dynamic Bioclimate Envelope Model (DBEM) with a size-based trophic model. The new approach provides spatially and temporally resolved predictions of changes in species' size, abundance and catch potential that account for the effects of ecological interactions. Predicted latitudinal shifts are, on average, reduced by 20% when species interactions are incorporated, compared to DBEM predictions, with pelagic species showing the greatest reductions. Goodness-of-fit of biomass data from fish stock assessments in the North Atlantic between 1991 and 2003 is improved slightly by including species interactions. The differences between predictions from the two models may be relatively modest because, at the North Atlantic basin scale, (i) predators and competitors may respond to climate change together; (ii) existing parameterization of the DBEM might implicitly incorporate trophic interactions; and/or (iii) trophic interactions might not be the main driver of responses to climate. Future analyses using ecologically explicit models and data will improve understanding of the effects of inter-specific interactions on responses to climate change, and better inform managers about plausible ecological and fishery consequences of a changing environment.

Comparison of new and primary production models using SeaWiFS data in contrasting hydrographic zones of the northern North Atlantic.

The accuracy of two satellite models of marine primary (PP) and new production (NP) were assessed against 14C and 15N uptake measurements taken during six research cruises in the northern North Atlantic. The wavelength resolving model (WRM) was more accurate than the Vertical General Production Model (VGPM) for computation of both PP and NP. Mean monthly satellite maps of PP and NP for both models were generated from 1997 to 2010 using SeaWiFS data for the Irminger basin and North Atlantic. Intra- and inter-annual variability of the two models was compared in six hydrographic zones. Both models exhibited similar spatio-temporal patterns: PP and NP increased from April to June and decreased by August. Higher values were associated with the East Greenland Current (EGC), Iceland Basin (ICB) and the Reykjanes Ridge (RKR) and lower values occurred in the Central Irminger Current (CIC), North Irminger Current (NIC) and Southern Irminger Current (SIC). The annual PP and NP over the SeaWiFS record was 258 and 82 gC m-2 yr-1 respectively for the VGPM and 190 and 41 gC m-2 yr-1 for the WRM. Average annual cumulative sum in the anomalies of NP for the VGPM were positively correlated with the North Atlantic Oscillation (NAO) in the EGC, CIC and SIC and negatively correlated with the multivariate ENSO index (MEI) in the ICB. By contrast, cumulative sum of the anomalies of NP for the WRM were significantly correlated with NAO only in the EGC and CIC. NP from both VGPM and WRM exhibited significant negative correlations with Arctic Oscillation (AO) in all hydrographic zones. The differences in estimates of PP and NP in these hydrographic zones arise principally from the parameterisation of the euphotic depth and the SST dependence of photo-physiological term in the VGPM, which has a greater sensitivity to variations in temperature than the WRM. In waters of 0 to 5C PP using the VGPM was 43% higher than WRM, from 5 to 10C the VGPM was 29% higher and from 10 to 15C the VGPM was 27% higher.